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In a recent letter we have presented the bare finite temperature effective po- 
tential for regularized SU(2) Yang-Mills theory taking the Haar measure into 
account H] . Here we present the renormalized effective potential at the order 
g 2 and show that it develops a non-perturbative minimum for sufficiently 
strong coupling, i.e. below a critical temperature. We discuss that this min- 
imum corresponds a phase which can be a candidate for the confining phase 
of the continuum theory. 

As to confinement in pure Yang-Mills theories general arguments were 
given that the deconfining phase transition occurs with increasing tempera- 
ture if the theory is confining at zero temperature || |3|]. It was also shown 
that lattice gauge theory does not confine static quarks if the Haar measure 
is replaced by the Euclidean one ||]. Non-trivial phases for SU (2) continuum 
theory has been predicted in Ref. |§. 
■ The 2-loop contribution of the order g 2 to the effective potential of SU (N) 

\q . gauge theory is well-known in the perturbative regime || [/]]. Now we deter- 

mine the effective potential for 577(2) up to this order by using the non- 
perturbative mean-field approach and treating the fluctuations around the 
mean field perturbatively. This mixture of the non-perturbative and per- 
turbative approaches provides some insight in the mechanism producing the 
Ph. non-perturbative phase below the critical temperature. We assume that the 

r-| ! component A^ =0 a=3 (// Lorentz index, a SU (2) colour index) of the gluon 

vector potential has the non-vanishing expectation value, v = (a'^^figA 03 ) (a 
the lattice spacing for the regularized path integral, /3 the inverse tempera- 
ture). Doing so we cannot assume any more that the field A 03 is small and, 
consequently, we must not replace the Haar measure, 
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(with a§ = gfia^Af, |0|, y]) by the Euclidean one in the path integral . 
This is in contrast to that one generally does in the perturbative approach. 
Lattice results have shown that confinement and global center symmetry 



are strongly related JK|. Since global center symmetry of the path integral 



is guaranteed by the usage of the Haar measure, the question arises whether 
does the Haar measure influence the effective potential. There are general 



1 



arguments that it does not influence the physics in the perturbative regime 



12, 13 1 . It was shown for SU{2) Yang-Mills theory that the tree level contri- 
bution of the Haar measure is cancelled by a piece of the 1-loop contribution 
of the longitudinal gluons . Thus the Haar measure does not play any role 
in the effective potential at the order g° even if a finite background field is 
assumed. In |J we have shown that the Haar measure reveals itself in the 
bare effective potential of the pure SU(2) Yang-Mills theory at the order g 2 . 
The authors of || found a cancellation of the Haar measure induced terms 
using a different approximation. Later we discuss the possible reason of this 
disagreement. Nevertheless, our results for the bare theory are in qualitative 
agreement with the findings for lattice gauge theory || as discussed in [p]]. 
Here we show that the basic features of the renormalized effective potential 
are not influenced by the Haar measure. It can reveal itself only in the fine 
details of the effective potential (in the approximation used). 

Our strategy of treating the IR and UV divergences is as follows. The 
bare theory is formulated in terms of three dimensional parameters: the 
inverse temperature j3, the UV momentum cut-off A, and the IR momentum 
cut-off fi. Thus UV divergences can occur as powers and logarithms of (3A. 
and A//i. UV divergences of the first type are those controlled by power 



counting and are removed by subtracting loops in the limit j3 — ■> oo JL4 
Of the second type UV divergences accompanying the IR divergences need 
some particular treatment. Their counterterms are introduced by choosing 
the IR momentum cut-off as a power series of 1/A and the logarithmically 
divergent terms are neglected. This removes both IR and UV divergences. 
The IR singularities are now removed in a way getting independent of the UV 
cut-off A when approaching the continuum limit (the rather correct approach 
instead of our earlier one [Q ) , since the IR problem in the real world is solved 
via dynamical generation of electric and magnetic masses, defining a physical 
mass scale independent of the UV cut-off. 

In M the effective potential was determined at the order g° allowing 
for the non- vanishing vacuum expectation value C/g = (a//3)(v/g). In the 
present paper we determine the terms of the order g 2 of the effective potential 



by using the general techniques of finite temperature field theory [[TJ]]. Fol- 
lowing U we use a time independent, diagonal gauge, and periodic boundary 
conditions for the spatial components of the vector potential in the 'time' di- 
rection. On the other hand we take all the terms of the order g 2 with in the 
effective potential. In this gauge the Haar measure induced potential in the 
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tree level action takes the form: 

1 ( 

Vntree = J d?X 111 (l - COS « 3 (f)) . (1.2) 

Our discussion starts with establishing the 2-loop corrections in terms of 
the gluon propagators. Then we determine their explicit expressions reveal- 
ing their dependence on the IR momentum cut-off /i and the UV momentum 
cut-off A, and perform the first step of the renormalization procedure by 
subtracting the corresponding loops at zero temperature. As to the next 
we complete the renormalization procedure by removing the additional UV 
divergences together with the IR divergences. Finally, we discuss the phase 
structure of the continuum SU(2) Yang-Mills theory on the base of the renor- 
malized effective potential obtained. 

Throughout this paper we use the notations of Ref. || and use the same 
cut-off regularization. The UV cut-off A is interpreted in terms of the lattice 
spacing a used for the definition of the path integral via a~ 3 = (2tt)~ 3 J a d 3 k 
with A = (67r 2 ) 1 / 3 a _1 and d 3 k the volume element in 3- momentum space. 
The same spacing a is assumed in the 'time' direction. 



2 The bare effective potential at 2-loop order 



Let us allow for the non-vanishing expectation value (a//3)(v/g) of the com- 
ponent A 03 of the gluon vector potential and shift the field variable by it. 
Inserting A 03 (x) = (a/ /3)(v/g) + 8<f>(x) in the Euclidean action 



S YM = \t ( dx° [ d 3 xF flu a F» v a (2.1) 

4a 4 n jo J 

of the Yang-Mills system (_F^ a the field strength tensor), the tree level action 
takes the following form: 

Stree = Sq + S\ + S*2. (2-2) 

The term 



2a 4 h Jo 
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dx Id x a z 



{V l 5<f)) 2 + (d°A ia + {v/(3)e a3b A ib 



+-(d i A ja -d j A ia ) 2 



[ d 3 xhi(l - cosv) + J_^o^_ f d 3 x ^2 / 2 _ 3 ) 

J 2a 6 1 — cosf J 
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contains the action of the free gluon field of the order (gy/h)° and the Haar 
measure induced self-interaction of timelike gluons, a term of order g 2 h (the 
fields are of order h 1 ^ 2 ). The other Haar measure induced vertices are ne- 
glected. Si and S 2 contain the usual 3-gluon and 4-gluon interaction vertices: 

Si = -L-Jd 4 xa [e aM 5<P A ld (d Q A ia + (v/ (3)e a3c A lc ) + e abc A ib A jc d l A ja } , 

(2-4) 



S 2 = JL^ J d 4 x ^ ca e Ma A ic A id {5(t)) 2 + l e ^c e ade A ib A jc A id A je^ ^.5) 

(e abc the completely antisymmetric tensor of rank 3). The self- interaction 
terms Si and S2 are of the order gh 1 ^ 2 and g 2 h, respectively. 

The main goal of this section is to determine the effective potential as 
function of the constant background field v including all the terms being of 
the same order g 2 as the one-loop contribution from the Haar measure. Then 
the interaction Si and S2 must be taken into account at 2-loop level. Our 
procedure is the following: 

• determine the partition function by treating the interaction terms Si 
and S 2 perturbatively; 

• determine the effective potential from the 1PI diagrams. 

In order to calculate the perturbative corrections we introduce the parti- 
tion function in the presence of the external sources j m and q: 

Z[f a ,q] = exp{-Si-S 2 }z [r,q] (2.6) 

with the generating functional of the free fields 

Z [j,q] = J V8<f>VA ia e- s °-^-^ =Zo[0,0]z[j,q], (2.7) 

where 

z\j,q] = expfyjDj) + ^(qVq)} (2.8) 
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and 



Z [0,0] = (DetM -Det(N + N H )y 1/2 exp^a- 3 j d 3 x\n(l -cost-)} . 

(2.9) 

Here D = M^ 1 and V = (N + Mh)^ 1 are the free propagators of the spatial 
and timelike gluons, respectively, defined through 

S = -a' 3 I d z x\Yi(l-cosv) + -(AMA) + -(5(j)(U + U H ) 5<P). 

(2.10) 

The matrices D and M. (V and J\f, Mh) possess Lorentz, colour and space- 
time coordinate (spatial coordinate) indices, (. . .) is a shorthand for contrac- 
tion. The operator S int — S± + S 2 is defined by replacing the fields A ia and 
5<p by functional derivatives 



fijia ' 
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(2.11) 



with respect to the corresponding external sources in the expressions for Si 
and S 2 - 

We expand the exponential operator in powers of (gh 1 ^ 2 ): 



-s 2 + \(s 1 ) 2 



+ 0((gVh 



e~ bint = 1-S 1 + 
After some algebra we obtain (h = 1): 

Siz\j, q] = -gjdxa [e 3da (Vq) [{Dj) id (d°(D 3 y a + (v/ ' (3)e 3ca {Dj) ic ) 



(2.12) 



3ca T^iidc 

xx 



+e 



+d» x DZ? + (v/(3)e 3ca D 

H)jy'v).,y'ir r n)jy" + Df x c d\D 3 y a 

HDjf^D"? + (Djy%D^ } z[j, q], (2.13) 



Siz\j,q] Q = 0, 



(2.14) 
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and 



dx |i 



2 9 



I ^abc^ade 



3ca 3<io«T^ 7-) 

C C i-^ J-S 



iicd 



^jjijbcjjijde _|_ jjijbejjjicd _|_ jjiibdjjjje 



}.(2.15) 



The 1PI part of S'i(5'iz)|o can easily be identified by noticing that (i) the 
terms containing first or third functional derivatives of z[j m , q] with respect 
of the external sources vanish, (ii) the non-vanishing terms with second func- 
tional derivatives of z[j m ,q] are 1-particle reducible. Thus the terms not 
containing the derivatives of z are the only ones contributing to the 1PI 
part: 



-Si(Siz) 



where 



olPI 



'VDD' + 'DDD', 



(2.16) 



'VDD' 



"ida Jieq 9 
£ £ 

2a 8 



d 4 yd 4 xV 



xy 



■a 



r^kiedpp pp r^kiga , aO y)kiea _ aO y\kiga 

xy x y xy ' y xy x xy 



+ {v/Py f9 (D^ d d° y D%f a + D%f d d° y D% a ) 
+(v/(3)e^ (D^ d dlD k ;f + D^dlDff) 
+ (v/P) 2 e 3ca e 3f9 (D^D k ^ c + D k xy ec D k xy fd 



(2.17) 



'DDD' 



2a 8 



e abc e gef 



d A yd A xa 2 



jjkieb f)ljfcf)i flk r)ljga 
xy xy u y u x xy 



1 pkieb^i j^ljfa _ ofc j^ljgc , j^kjec plifbfli ofc j^ljga 

' xy y xy rr. rr.il ' rr.it rr.it ^ll^rr. rr.it 



xy 



xy 



xy 



y^x xy 



1 r\kjecryi j^ljfa _ ofc j^ligb , j^lifbfli r^kjea^k f)ljgc 

' xy u y xy u x xy ' xy u y xy w x xy 

1 Tfii c (f J~) k j ea f) k f) li 9 

' Xy y xy ^rr. rr.it 



■HI 



(2.18) 



The effective potential is defined by the 1PI part of the partition function. 
The 1-loop contribution of the gluon self-interaction Si nt to the effective 
potential vanishes identically due to Eq. ( |2.14j ). The 2-loop contribution 
AV e ff is given by 



Z[0,0] 1PI = Z [0,0]exp{-/3VAV e// } 



(2.19) 
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Figure 1: 2-loop contributions to the effective potential. 



with 



AV eff = (m- 1 ^43,Q}\ -^(Si) 2 43,Q] 



olPI 



h + h + h + h 



(2.20) 



(V — > oo the 3- volume of the system). The various terms on the r.h.s. of Eq. 
( |2.20| ) correspond to the diagrams in Fig. 1 and represent the bare 2-loop 
contributions to the effective potential: 



./': 



IB 



3_ 

2a 4 



e abc e ade / )'J'J' + Df^D^f + , (2.21) 



= V-V (D 

u xx 



iicd 
xx 



y ;:>:-> 



D l , 

XX XX I 1 



(2.22) 



3B 



1 

2(3Va s 



g 2 e abc e gef a 2 j \ D^f D^d^D^f 
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ijjkiebfli Tjljfa . aky)ljgc , y)kjec jjlifb oi ofc Wjga 

' xy x xy x xy ' xy xy w x w x xy 

I jjkjecfli Tjljfa . nfc r^ig6 

~ xy x xy x xy 

i JjUfbfji T~)kjea _ afc r^ij'gc , jjljfcpji f\kjea _ afc rjiigfe 

' xy x xy x xy ' xy x xy x xy 



(2.23) 



IB 



i. e 3da e 3eg f ^4,^4 



2f3Va 



d yd xT> xy ■ 

r^kied / pp\2 -r^kiga , p,kiea / cp\2 p>kigd 

xy v^x) xy ' xy \^x) xy 

+ {v/Py f9 {-Dxy ed d° x D k xy f a + D% a d°D% d 
+ {v/(3) 2 e 3ca e 3f9 (D k xy ed D k xy fc + D k xy ec D k xy fd 



(2.24) 



Later we shall discuss how to obtain the renormalized contributions f\, f'2, 
fs, and fi from the corresponding bare ones. 

The 2-loop corrections ( |2.21| )-( ]2.24| ) are given in terms of the free gluon 
propagators. After changing to momentum representation we read them off 
from the free action Sq. Then we obtain for the free propagator of timelike 
gluons 



V(k) 



a 2 (k 2 + M 2 



(2.25) 



with M 2 a 2 = a g [3(1 — cost>) 1 and for that of the spatial gluons 

D ijab (u n ,p) = a~ 2 \(u 2 + p 2 )5 ij 5 ab + {v/(3) 2 5 ij {5 ab - 5 a3 5 b3 ) 



- P Y6 ab - 2(v/(3)iu n e 3ab S lJ 



(2.26) 



with the bosonic Matsubara frequencies u n = 27in/j3 [n = 0, ±1,±2...). 
Separating the transverse and the longitudinal parts the matrix on the r.h.s. 
of Eq. ( |2.26p can be inverted: 



D ijab (u n ,p) 



pipj 



a- 2 A ab (Lu n ,0)^ + a- 2 A ab (u; n ,p} I 8» - 



pipj 
p 2 



(2.27) 



S 



with 

/ imp)+d-(p)) -\{d+{p)-d-{p)) o \ 

A a& K,p) = \+Wm-d-{p)) \{d+{p) + d-{p)) 

V d n (p)J 

(2.28) 

and dt(p) = l(ut? + V 2 }-\ d n (p) = (u 2 n + p 2 )-\ ^ = uj n ± (v/(3). In 
coordinate representation, the free propagators are then given by 

V 3 r d 3 k exp{ik(x-y)} = a exp{-M|f - y\} 

Sff " M27T) 3 a 2^2 + a 2 M 2 "27T |f-j7| ' ^ ' ^ 



D%f = /3 _1 E/ -^D^^p)^- 3 '^' 1 ^-^. (2.30) 

In order to avoid IR divergence we introduce the IR momentum cut-off fi 
and replace A ab (u n ,0) by A ab (u n ,fi). 



3 Explicit form of the 2— loop contributions 

The explicit expressions for the various contributions f a (a — 1, 2, 3, 4) to the 
effective potential were obtained by (i) performing the contraction of colour 
and spatial Lorentz indices; (ii) performing the Matsubara sums by contour 
integral technique; (in) expressing the result in terms of the gluon occupation 
numbers. The first step of renormalization was carried out by subtracting 
from each Feynman diagram the similar ones with a loop taken in the limit 
(3 — > oo in all possible ways, and neglecting all terms of the effective potential 
independent of the temperature and of the background field v. Below we shall 
see that this renormalization procedure removes the terms which are linear 
in the gluon occupation numbers. Exceptions are some terms containing the 
occupation numbers of the zero mode (that for gluons with the IR cut-off 
momentum //). Therefore the 'renormalized' contributions f a obtained in this 
chapter still contain the UV divergences accompanying the IR divergences. 

We analize the contributions f a (a = 1,2,3,4) and separate the pieces 
f a g i and f a vac corresponding to the gluon gas and the vacuum, respectively: 
fa = fagi + fa vac- Each of these terms can be written as the sum f a ... = 
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Figure 2: 'Renormalized' contribution fi to the effective potential. 



fa... + fa... where /° and are the contributions for v = and the 
additional contributions due to v ^ 0, respectively. The terms depending 
on the occupation number of the zero mode are considered as part of the 
vacuum contribution. 

Below we do not write out the terms of the order A s with s = 0,1,2,3 for 
each contribution f a separately. However, we shall include all these terms in 
our final expression for the bare effective potential in Sect. 5. 

• Contribution fx 

The 'renormalized' contribution fi corresponds to the diagrams in Fig. 
2. Performing the subtractions on the r.h.s. of Eq. ( |2.21| ) we obtain 
for the 'renormalized' contribution 

h = -^<f <•""<'"'■ (lyj'riyj:! 1 ■ • »t) (3.1) 

with 

D i j x b = I)'?:' 1 ' - 2 lim Di Jx b - (3-2) 



The calculation described in Appendix A gives the following expression 
for the renormalized glue gas contribution for v — (see Eq. ( |A.14| ) ): 

fO _ , 2 [ d 3 pd 3 q n p n q _ g 2 

h ° l ~ 49 J (2tt)6 pq "36^ ^) 
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with the occupation number n p = [exp([3p) — 1] 1 of the gluon state 
with momentum p. Eq. ( [A. 15 ), leads to the vacuum contribution for 
v = 0: 



a 



1 a 2 



a 2 n 



167T 4 fid 47T 2 

1 



+ 



(2n„ + l)(2n. 



(3.4) 



3/3 2 /ia 4/z 2 a 2 

From Eq. (|A.16|) we find for the additional contributions due to v ^ 0: 



f A 
Jl gl 



3vr 2 /3 4 



V 




V 


2tt\ 


sin — 


( 


sm — 




2 




2 


~ ~3/ 



(3.5) 



J 1 vac 



9_ 
3 



2 sin 



+ 



Tra 2 (3 2 (fia) 
1 



(iV» + rv) + 



a 4 (/ia) : 



3a 2 /? 2 (/xa) v / ' 
(JV» - n M ) (iV» + 3n M ) 



(W» - n, 



(3.6) 



in terms of the 'generalized occupation number' N^v) = N p=fl (v) de- 
fined via Eq. (|A.7| ). The terms fi gl and f^ gl contribute to the energy- 
density of the blackbody radiation, whereas the vacuum contributions 
fi vac and fi vac are IR divergent and depend on the temperature. Their 
dependence on the IR momentum cut-off \x is revealed using the expan- 
sions: 



n u ps 1 8a \8aY -\ 

^ 8n 2 YT^ 720 VPW 720-42 



m 5 + 0(S)(3.7) 
(3.8) 



with A = (1 — cos-u) 1 , B = A — 3A 2 . Then we obtain: 

2g 2 A 6 



f° 4- f A 

J 1 vac 1 J 1 



(2tt) 4 27 

8 n 



2 2 , 
+ ^4 



»y 6 + 2A ) + ^[l A - A * 



+ C(A 3 )(3.9) 
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Figure 3: 'Renormalized' contribution f 2 to the effective potential. 



after removing the constants independent of the temperature and the 
background field v. Thus only the terms quadratic in the occupation 
numbers of the zero modes contribute. 

Contribution f 2 

Making use of the expressions (|2.25|) -( f2.30|) of the free propagators, we 
rewrite the bare contribution Q2.22|) as follows: 

,2 



./': 



g 



2B 



-Ma 



2vra 4 



a(fi,v) + 2R(v)\ 



(3.10) 



where a(p, v) is given by Eq. ( |A.6| ) and from Eqs. ( |A.9| ), and ( |A.10| ) 

we get: 



R(v) 



d 3 p 



a(p,v) 



+ 



a 



12(3 2 8vr 2 2vr/5 2 



v 

sin - 

2 



(3.11) 



The 'renormalized' contribution is represented by the diagrams in Fig. 
3 and is given by 



9 _ e ~Ma 



2vra 4 



9 _ £ -Ma 



2vra 4 



a 1 (a(n, v) - (Tooi/i, v)) + 2 (R(v) - Roo{v)) 
jj,a 



a 2 




i 


V 






Q- 




sm - 


)] 


IP 




7T 


2 





(3.12) 



with Ooo(yu, v) = lim ^oo cr(/ ^, v) = OoJjj,, 0) , Roo(v) = limbec R(v) = 
Roo(0) given by Eqs. ( |A.13| ) and ( |A.12| ). Here we made use of the fact 
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that the loop with the propagator T> vanishes in the limit (3 — > oo due 
to M ~ (3 — > oo. The contribution / 2 is IR divergent and depends on 
the temperature. It can be considered as part of the free energy density 
of the vacuum, f 2 vac = /2 an d f% g i = 0. For v — 0, i.e. M — > oo it 
vanishes, /fmc = 0> an d f° r " ^ 0, i-e. finite M we get /^ac = fa- 
Using the expansion ( |3.8| ) we obtain: 

ftac =h = -^L\^p a + ( 3 - 13 ) 



Contribution f% 

After Fourier transformation the bare contribution f 3 s given by Eq. 
(2.23|) takes the form: 



?>B 



2a s 



g2 e abc e gef ft-3 



d 3 pd 3 qd 3 Q 
(2tt) 9 



■(2tt) 3 5(p + q + Q)(35 ni+m+na ■ 
■a 2 ^D kieb (u ni , q)D ljfc (u n2 , p) 

+D*"(u ni ,$D lifb (u na ,d) Q i Q k D^\u n . i) Q) 

+ (D^\u ni ,q)D^ a (u; n2 ,p) 

+D u f\u ni ,q)D k i° a {u n2 ,if)) p i Q k D lj9C (uj n3 ,Q) 

+ (D k ^ c (uj ni ,q}D l ^ a (uj n2 ,p} 

+^>(o; nu ^D^ eo (a; n2 ,^)p i Q fe ^ 6 (a; n 3,g)] • (3.14) 

The 'renormalized' contribution of the diagrams in Fig. 4 is given by 
(Appendix B): 



fs 



-2 9 



d 3 pd 3 Q 
(2vr) 6 



1. 



where the integrand takes the form: 

I = 1 2 + T 2 + l x + T[ + Ip + Tp + J . 



(3.15) 



(3.16) 
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Figure 4: 'Renormalized' contribution fa to the effective potential 



— # 

With the help of the momentum dependent functions P<p s (q,p,Q), 
$>(q,p,Q), and 4> (q,p,Q) defined in Appendix B by Eqs. ( P.14j) , 
( |B.15| ), and ( |B.9| ), respectively, and the 'occupation numbers' = 
[exp(/3£f ) - l}- 1 with Ef = q±iv/fl we find the expressions below for 
the various terms of the integrand. (For equalities holding only up to 
the terms containing odd powers of the scalar product pQ stands ~.) 
The terms quadratic in the occupation numbers give contributions to 
the blackbody radiation term: 



1, 



-{n; + n p 



-(n+ + n_) 



+ n^ + n Q 
4Qp 

4wq + + Uq 
4Qp 

4n Q + Uq + rig 4 
4Qp 3 



1 (pQ)Pfo P 2 + Q 2 
2(pQ) 2 -p 2 Q 2 (p 2 -Q 2 ) 2 



1 + 



p 2 Q 2 



(3.17) 



T 



-g(n+-n p )(nQ--n Q ) — 



(pQ) 2 - P 2 Q 2 



(3.18) 



The term X\ is linear in the occupation number at the IR cut-off and 
is IR divergent: 



Ji 



2n Q + nZ + n Q n+ + n 



2Q(p 2 - Q 2 ) 



+ 



n Q+ n Q n n 

Q(p 2 - Q 2 ) fi 
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P 3 II 



[n+n p + 2n p ) (raj + ) + 2(n+ + n p )ra J o (g, P, <?) 



(2n Q + + ra Q ) (raj + ra^ ) + (nj + 



Q 3 2Q(p* - g 2 ) 



(3.19) 



The term X( vanishes for vanishing background field identically; it is 
IR finite: 



T , _ 3 n+-n Q 



2 (p 2 - Q 2 ) 2 
2 



n 



t-n-)mviQ) 



n 



Hq,P,Q) 2 ^ 

(p2_Q2 )2 + Q4^(g.P.Q) 



(3.20) 



The terms quadratic in the occupation number at the IR cut-off are 
given by: 



1 



2p 2 fi 2 



[n+ + n~ + 4n„)(raJ + n~)MQ,P, Q), (3-21) 



and 



= "^("J -^J 2 0o(?,P,<2)- 



(3.22) 



The integrand X/3 is temperature dependent and IR divergent. Finally, 
there is the IR divergent, temperature independent part of the inte- 
grand: 



Xn = 



1 



2qpQ(p + q + Q) 



3<f>(q,p,Q) (I 1 



p 2 jj? 



3 ^-^W,p,q). 
p j 



2pQ(p + Q) V/i P + <2/ 

(3.23) 
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The corresponding contribution to the effective potential is a part of 
the energy of the perturbative vacuum. 

Introducing the notation 
we write 

h = hi + + hi + fk + fw + hi + ho- (3.25) 



Here the terms fa, fa, and contribute only for u^O, whereas the 
other terms can be cast into a piece for v = and the rest vanishing 
for v = 0: fa = / 3 °2 + /si /si = / 3 °i + /£, and = f% + f 3 % 

The only terms not suffering from IR divergences are fa and fa, the 
contributions to the density of the free energy of the glue gas, 



h 9 i = & + f£ + fL- (3-26) 

For v = we find: 



hgl-J32 ~ ^9 / /„ 3g 77" - 79/ o 4 - 



d 3 pd 3 Q nQn p g 2 
(2tt) 6 ~ 72> 

It is a contribution to the energy of the blackbody radiation. 

For v 7^ we obtain the additional contribution to the density of the 
free e nergy of the g lue g as, f^ gl = f& + fa where (see Eqs. ( |B.36| ), 
(|R37]), and (|R38l) -( p^0l) in Appendix B) 



h 



7 



sin 2 v 



32 



2e- 2 (2vr) 4 f3 A (cosh 1 - cosv) 



(l-(l/e))(2-(5/e)) 



(3.29) 



with G(v) defined by Eq. (|B~37|) . 
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The IR divergent terms 



h vac — /31 + hp + /30 + /3I + f'z 



3/3 



(3.30) 



contribute to the free energy of the vacuum. Their dependence on the 
IR cut-off /i is revealed as described at the end of Appendix B. Thus 
we obtained ( A = (1 — cost*) -1 ): 



2<? 5 



1 A 7 25 



3 vac 



(2tt) 4 \/i 4 (3 128-7 
93 



/i 3 



1 

? 
1 

+— 

/i 2 

+A 7 /3 
-A 6 



„ A 6 2431 A 6 , A 

A H 1 In — 

(3 16 • 9 • 5 • 7 12/5 /i 



-A 7 - 

128-3-7 

A^ 33 / o \ A 6 4, A 
732-^ ln( ^-T2 A1 ^ 
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16-9-7 
10039 
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32e 
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64 



64-9-5-7 



A 



A 5 
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3 e 



+■ 



(A - 3A 2 ) - A 5 



64 -3 V 7 64-5 

1 A 6 /5 /l \ A 
[ - + 2A ] In - 

A* 



ln(/^) 



1 

+- 
A* 



/i 24 V6 ' "7 
A 7 /? 2 /5 _ 4 
3~e 



+A 6 /3 
+ 



64 
2431 



32-9-5-7 V6 

A 7 /? 3 / 1 



93 



+ 2A 



281 
4-3-7' 



a 5 1423 A 
A 5 A 

64-3-5 



3A 4 



4/3 



+ 



-- + A 2 ) - 
360 3 

A^ 



64 
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48 V ; 16 



A 5 /3 /13 1531^ 
~30~ VT + ~32~ 

In* 



ln(/3/x) 



A 7 /? 3 
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Figure 5: 'Renormalized' contribution f± to the effective potential 



+ f ™» - 5807.4 

128 ■ 9 ■ 5 ■ 7 V 2 

A „„ / 4851 159 1587 „ 

+A 5 /? 1 1 A 

H V 128-9-5 16 • 3 • 5 e 128-9 

77 . 23 



-A + —A 



32 • 3 e 24 
551 

+A 4 4 o ' - \ + G(A 3 ). (3.31) 



4-3-5-7. 

The terms independent of (3 and v are neglected. 
Contribution f± 

As discussed in Appendix C we obtain the expression ( |C12| ) for the 
contribution f± after carrying out the subtractions shown in Fig. 5. 

In the perturbative regime with v = 0, i.e. M — > oo we get / 4 = f% mc = 
from Eq. ( C.12 ). This is because the infinitely heavy timelike gluons 
do not propagate. 

Following the procedure described at the end of Appendix C, we get 

= U = "(SiRiS^J + 2 f) + ° (A3) - (M2) 

Here the terms independent of (3 and v are neglected. 
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4 Blackbody radiation 



The complete 2-loop contribution to the density of the free energy of the 
glue gas is given by f gl = /°, + 

Consider at first f gl for the case v = 0. The renormalized diagrams 2 and 
4 (Fig. 1) do not contribute to the glue gas as discussed in Sect. 3. The 
renormalized diagrams 1 and 3 (Fig. 1) give the gauge invariant contribution 
to the free energy density of the glue gas (see Eqs. ( |3.3|) and (|3.27| ) ): 



f& = /i% + / 3 % = (4 + % 2 / 



d 3 pd 3 Q nQfip g 1 



(2tt) 6 Qp 24/3 4 



(4.1) 



This result is just the same as that one obtains by adding the contribution of 
the gauge field and that of the ghosts in the usual perturbative calculation 
given in M: 



ck9 



d 3 pd 3 Q tiq 



(2vr) 6 Qp 



(4.2) 



with ck = c g host + ci + C3, where the weights of the ghosts, diagram 1 and 
diagram 3 are c g h os t = c\ = -f'6, and C3 = — |-6, respectively, for SU(2), 
so that ck = Q- The contributions of the single diagrams are different in our 
case, they are gauge dependent, but the sum of them is gauge invariant, as 
expected p. 

For v 7^ the additional contributions to the free energy of the glue gas 
are due to diagrams 1 and 3, as f^ gl = f£ gl = 0. According to Eqs. ( |3.5|) , 
O, and flpD we find: 



Jgl 



fl gl + ^3 gl 

24/3 4 



" 8 


V 




V 


2tt\ 




sin - 


( 


sin - 




7T 2 


2 




2 


~ ~3J 



G(v) G(v) 



21 



sin 2 v 



47r 4 e 2 (cosh 1 — cosv) 




;i-e- 1 )(2-5e- 1 ) 



(4.3) 



with e the basis of the natural logarithm, and the function G(v) given by Eq. 
(p.37 ). This contribution vanishes for v — 0. 

The sum f g i = f gl + f£ represents the complete contribution of the glue 
gas to the free energy density. This is an IR finite contribution, independent 



19 



of the IR cut-off momentum //. The free energy of the vacuum, however, 
depends on the background field v and the IR cut-off momentum, as well. 
In order to reveal its v dependence, we have to get rid of the dependence on 
the IR momentum cut-off. 



5 IR divergent contributions 



Gathering the IR divergent pieces obtained in Sect. 3 we find the following 
contribution to the effective potential: 



V IR 

~ fl vac + fl vac + fl vac + f^ vac + f, 

= V 2 + V 



A 

4 vac 



with A = (I — cost") 1 and 

2g 2 



(5.1) 



+/i" 3 [-A 7 b 31 A + A 6 f3-% 2 - A 5 f3-\ 3 v' 
A 6 

+/i ~ 3 l2^ ln (A//U) 

+/i" 2 [A 7 (3(b 21 - b 22 A - b 23 A 2 ) - A% 24 A 

+A 5 /T 1 (-& 25 + b 26 A) - A A (3-% 21 v 2 
+fi~ 2 [A 5 (3-% 8 hitfn) - A 6 (A/I2) ln(A//x) 

+A 3 p- 3 T 23 (v) + A[3- 5 T 21 (v) + r 6 r 2Q (v) 
[A 7 f3 2 {-b ll + b l2 A + b l3 A 2 ) 

+A 6 f3(b l4 + b 15 A) - A 5 b w A - A 4 (3-% 7 

-A 3 p- 2 T 13 (v) - Ap~ 4 T u (v) - (3~ 5 T w (v) 

[A 7 (3 2 b 18 (A - 3A 2 ) - A 5 (6 19 + h 10 A)] ln(/fy) 
[A 6 (/?/144) (I + 12A) + A% U A] ln(A//x) 

+A 7 p 3 (b 01 + b 02 A - b 03 A 2 - b m A 3 ) + A 6 f3 2 (-b 05 A + b 06 A 2 ) 

+A 5 (3(-b 07 + b 08 A + b 09 A 2 ) + A% W A 
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+A 3 /r 1 r 03 ^) + A 2 /r 2 r 02 (*;) + A(3- 3 r i(v) + /r 4 r, 



+ 



+ 



A 7 /? 3 /J_ 
64 V360 



A + A 2 +A 5 /?(6 011 + 6 012 A) 



00 [V ) 

A 3 

4/9 



ln(/3/i) 



A 6 /3 2 
48 
2.54 



(A - A 2 ) + A% 13 A - 0.52A 2 /T 



V 



(1 - 2A) 



ln(A//i) 



(5.2) 



and 



V'o 



i - I ^l 3 ^{^A3 r V^(MA/,)-1.4) 

+ A T 1 A 2 /rV 24tt 4 + 7T 4 (-2.25 A 3 + 78A/TV) ln(A//i) 
-643A 3 + A/r 2 (-130 + 1022t; 2 )} (5.3) 
where the constants 6 . are positive numbers: 



641 


= 0.028, 631 = 


0.035, 632 = 


= 0.408, 633 = 


0.086, 






hi 


= 0.0054, 622 = 


= 0.011, 6 23 


= 0.016, 6 2 4 


= 0.424, 625 


= 0.287, 6 26 = 


= 0.206, 


hi 


= 0.171, b 28 = 


0.206, 










bu 


= 0.003, 612 = 


= 0.023, 613 


= 0.052, 614 


= 0.033, 615 


= 0.408, &16 


= 1.48, 


617 


= 0.793, 6 18 = 


0.0052, &ig 


= 1.65, 61 10 


= 0.667, 61 11 


= 0.063, 




601 


= 0.0003, 602 = 


= 0.008, 603 


= 0.076, &04 = 


= 0.0078, 6 05 


= 0.108, &06 


= 0.107, 


&07 


= 0.60, 6 08 = 


0.66, 609 = 


0.96, 6010 = 


1.31, 6 ii = 


= 0.145, 6012 


= 1.59, 



60 13 = 0.063, 
and 



r 2 3^) 


= 14.35 - 0.23A + 0.054 2 - 2.79 


. v 
sm — 
2 


(5.4) 


F 2 i(v) 


= -3.48 + 0.284 - 0.0544 2 , 
= -0.60 + 0.0484 -0.0164 2 , 




(5.5) 
(5.6) 


ris(u) 


= 17.66 - 0.234 + 0.054 2 - 2.79 


. v 
sm - 

2 


(5.7) 


rn(u) 
rio(u) 
r 03 (^) 


= -0.52 + 0.284 - 0.0544 2 , 
= -3.24 + 0.164 -0.0164 2 , 
= 1.429 - 3.964 -0.1824 2 + 0.0254 3 


(5.8) 
(5.9) 



-(0.233 + 1.404) 



. v 
sm - 
2 



(5.10) 
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T 02 {v) = 1.94 v 2 , 

r 01 (v) = -0.064 - 0.904A + 0.023A 2 



0.008A 3 . 



(5.11) 
(5.12) 



The contribution Vq of the order g° arises from the diagram f'4 through 
expansion in powers of (fi/M) 2 ~ /i 2 /3 _1 A~ 3 g~ 2 . Thus terms of the order 
g 2 (l/g 2 ) = 1 occur. The terms of higher order in (l/<? 2 ) vanish as some 
powers of (1/A). 

The bare finite temperature effective potential V e ff{f3,v) as function of 
the constant background field A 03 = (a/p)(v/g), A 01 = A 02 = A ia = 
(i = 1, 2, 3; a = 1, 2, 3) takes now the following form: 



where 



Veff(P,v) 



V w {/3, v) + V H (J3, v) + AV(P, v) + Vj 



IR, 



2if_ 
> 



is the term obtained by Weiss 



V, 



n 



1 

24 



ao_9_ 
a 4 sin 2 I 



I v I 

7T 



- 1 



(5.13) 



(5.14) 



(5.15) 



with cto = (67r 2 ) 1//2 is the term induced by the Haar measure, the 2-loop 
contribution due to the usual non-abelian self-interaction is given by Eq. 

as 



Jgl 



24/? 4 





V 




V 


sm - 


( 


sm - 




2 




2 



2tt 

T 



7T 



21 



sin 2 v 



47r 4 e 2 (cosh 1 — cosv) 



1 



-1 



3 

5e" x 



G(v)j G{v) 
(5.16) 



and the IR divergent piece is given by Eq. ( ^.1|) . As the variable A 03 is 
compact, we have to set forth periodically the expression of the effective 
potential given for v G [0, 2tt] . 

The terms periodic in v arose due to the longitudinal (spatial) gluons. 
As observed in P| their Matsubara frequencies are shifted by the background 
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field, u n ±(v Having rewritten the Matsubara sums over the longitudinal 
modes as contour integrals on the complex energy plane, they are dominated 
by the complex poles at \p\ ± i(v//3). Thus the generalized occupation num- 
bers rip and N p (v) arise and lead to terms periodic in v. 

Here we recover our observation made in [|IJ, that the Haar measure term 
~ +g 2 K A A of the bare effective potential is not cancelled by the appropriate 
IR finite piece ~ +g 2 b iqA 4 A due to the longitudinal gluons as 6 10 > 0. Our 
finding that the Haar measure term, g 2 A 4 A of the bare effective potential is 
not cancelled by a corresponding piece of the contributions of longitudinal 
gluons is in disagreement with the results of || . This disagreement can be the 
consequence of the difference of the both approaches. We treat the quadratic 
piece of the Haar measure induced potential (the last term of Eq. fl2.3|)) 
non-perturbatively, including it into the free propagator of the field <50, and 
the first term of Eq. (|2.5|) perturbatively. On the other hand both terms are 



treated non-perturbatively in ||. Below we show, however, that the Haar 
measure does not influence the main qualitative features of the renormalized 
effective potential in our approximation. Thus the question of cancellation of 
the Haar measure is not crucial as far as our conclusions on the basic features 
of the effective potential are considered. 

The bare effective potential represents a sum of IR singular terms, each of 
them being a sum of powers of the UV cut-off. IR divergences appearing as 
powers of the dimensionless parameter A//i lead to additional UV divergences 
of higher order than one would expect by power counting. These are yet 
present although the subtractions of loop integrals at zero temperature have 
been performed. 

Now we remove the IR singularities. For this purpose we choose the 
IR momentum cut-off fi as a power series of the inverse of the UV cut-off 
(i.e. of the lattice spacing). Actually we introduce the power series for the 
dimensionless IR cut-off in terms of l/(/3A): 

(Pfi)- 1 = X>(/3A)- fc . (5.17) 

k=0 

The series can be terminated by k = 7, as the terms of higher order will 
vanish in the continuum limit when the series is inserted in the bare effective 
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potential. For later use we write: 



{Pp)~ s = i> sfc (/3A)- fc , (a = 1,2, 3, 4) (5.18) 

k=0 



where the coefficients are obtained by taking the powers of (|5.17 ). The IR 



cut-off defined by Eq. ( |5.17| ) becomes independent of the UV cut-off for large 
values of A and proportional to the temperature T = 1/(3, p — > (l/p )T. 

For the sake of simplicity we neglect the logarithmically divergent terms. 
Then the coefficients on the r.h.s. of Eq. ( |5.17|) are defined by requiring 
that the coefficients of A 7 , A 6 , . . ., A vanish in the effective potential ( p-l^Q 
having inserted the sum ( |5.17|) in it. These renormalization conditions lead 



to highly non-linear equations for the determination of the coefficients p k 
(k — 0, 1, ... , 7). Assuming that our procedure gives a reasonable result, i.e. 
the inequality po 1 holds, the coefficients m s k in Eq. (|5.18|) take the rather 
simple form: m s o = Pq, m s k = sp^~ p, k f° r k = 1,2, ... ,7, and s = 2,3, 4. 
Keeping only the terms of leading order in p , the renormalization conditions 
give rise to the following equations linearized in p k (k — 1, 2, . . . , 6): 

= b 41 m 40 -b 31 Am 3 o + 0(nl), (5.19) 
= b 41 m 41 + b 32 m 30 + O{nl), (5.20) 
= b 41 m 4 2-v 2 b 33 m 30 + O(fil), (5.21) 

= &4im 43 + 0(/io), ( 5 ' 22 ) 

= 641^44 + 0(^0), (5-23) 

= h imi5 + 0(fxl), (5.24) 

= b A1 m 46 + 0(fil). (5.25) 

The approximate solution of these equations is given by 

// = PoA, pi = -pi, p 2 = v 2 p 2 , 
p k ~O(l/p )&0 for /c = 3,4, ... ,6 (5.26) 
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with p = &31/&41, Pi = ^32/(4641), and p 2 = 633/(4641). Thus the coefficients 
A*fc (k > 3) are suppressed and can be set zero in our approximation. As the 
series terminates at k = 3 it seems to be reasonable to choose fij = 0. 

In the continuum limit the IR momentum cut-off takes [l = for t> = , 
whereas it has its maximum value /i = 1.6/(3 ~ (tt/2)T for t> = tt. Thus the 
IR momentum cut-off is much less than the lowest non-vanishing Matsubara 
frequency 2-rT . This seems to be a justification for considering the coefficient 
/i large. 

Inserting the limiting value of the IR momentum cut-off in the expression 
( |5.13| ), we obtain for the renormalized effective potential: 



with 



V R = V W + V 0R + V 2R (5.27) 



V R = 0.097 w 4 /r 4 (l - cos v), (5.28) 



V 2R = ^ I ^{Po(v) + P 1 (v)A + P 2 (v)A 2 -OAOA i - 0.025A 4 

071 p t 

+ (-55.5 + 15. 7w 2 )| sin ^| +26.3 sin 2 ^j (5.29) 

where 

P {v) = 25.4 - 80.5v 2 + 0A6v 6 , (5.30) 

P x [v) = 35.5 + 1.29w 2 , (5.31) 

P 2 (v) = -3.9-0.28i; 2 + 1.27sin 2 w. (5.32) 

Here we used approximately G(v) ~ —0.32A — 1.2. The purely polynomic 
part of the renormalized effective potential arises due to the decay of the 
A 03 condensate into two transverse gluons and a timelike fluctuation 5$ 
which then annihilate reproducing the condensate. This is a piece of the 
contribution / 4 . The term ~ gv(A il A il + A i2 A i2 )5(j) in Eq. (£Oj) for the cubic 
self-interaction is responsible for this subprocess. 

Now we can come back to the question of the significance of the Haar 
measure as far as the renormalized effective potential is concerned. The 
Haar measure term Vh present in the bare effective potential explicitly as 
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found in [|IJ gives a contribution of the order one to the coefficient of A 4 , i.e. 
to the r.h.s of Eq. (|5.22 ). Thus it means a correction of the order 1///q to 
the leading term of /z 3 being of the order l//io- Thus we conclude that the 
Haar measure term does not influence the basic features of the renormalized 
effective potential in our approximation. Nevertheless it reveals itself in the 
details explicitly. 

The periodic terms of the renormalized effective potential occur due to 
the longitudinal gluons in our approximation, as discussed above. They are 
invariant under the global center transformation v — > (27r — v). The polyno- 
mials break the global center symmetry explicitly. This explicit symmetry 
breaking is the consequence of the ansatz (A 03 ) ~ v inserted in the tree ac- 
tion. It would be a way to preserve the global center symmetry that one 



makes use of the Yang-Mills action derived in |10 . 

Here we observed renormalizability at the order g 2 . Although the Haar 
measure generates infinitely many vertices, the general results for the non- 



linear sigma model |12j and those for renormalizable lattice field theories |fL3 
give the hope that the theory is renormalizable. 

At the end of this section let us make some remarks on the validity of the 
expression obtained for the renormalized effective potential: (z) V2 tends to 
infinity for v — 0, this must not be taken, however, seriously as the expansions 
used for the generalized occupation number are only valid for v 7^ 0. It means 
that our effective potential is valid in the middle of the interval v G [0,27r]. 
At the ends of the interval we have to accept the perturbative results of 
[0- B- Our approach of linearizing the renormalization conditions and 
solving them for the coefficients [ik of the expansion of the IR cut-off is 
rather crude as the coefficient /io becomes of the order one for v ~ 71". (Hi) 
It is unclear what effect the 'less dangerous' logarithmic divergences on the 
effective potential have. 



6 Phase structure 

Below we show that the renormalized effective potential exhibits a very re- 
markable feature. Namely, it has a minimum at non- vanishing background 
field v 7^ which becomes deeper than the minimum of V\y at v — describ- 
ing the perturbative phase. Furthermore, the position of the new minimum 
is rather close to v = it which corresponds to the vanishing Polyakov line 
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cos(u/2) = 0, i.e. the confinement. 

It is easy to see that the polynom Po(^) has a minimum for v pa 0.9n 
and the periodic functions in V^r have their minima at v = tt. Therefore we 
seek the extremum of the effective potential at v = tt — s with e <C tt. The 
necessary condition for having a minimum at v = tt — e is given by: 








dV w dV 0R dV 2R 



dv 



dv 



13.1e) + f-j (342.4 - 6534.3e) 



57P 



The solution is given by 



0.052 (g 2 + 53.0) (g 2 + 1.5 



(6.1) 



(6.2) 



leading to e = 1.8, 1.1, and 0.05 for g 2 = 0, 1, and g 2 — > oo, respectively. 
For strong coupling the position of the extremum is defined by the term V2R. 
For weak coupling the term Vqr defines the position of the extremum in our 
approximation. However, then tt — v gets large and our expression for Vr 
is not valid any more. It is easy to show (with accuracy of 0(g 2 ) ) that we 
obtain a minimum at v 7^ for any coupling, as 



d 2 V R 



dv 2 



/r 4 (52.7 + 9.2/) > 0. 



(6.3) 



It is well-known that there is a minimum at v = 0, since for weak coupling, 
i.e. small values of v the formula Vw modified by the 2-loop correction 
(1.1) is valid ||. Our expression is not valid for small values of v. For suffi- 
ciently strong coupling the minimum find by us is in the validity range of our 
expression for the effective potential. Then there should be a maximum some- 
where between the minima at v — and v = n — e. The value of the effective 
potential at v = is given by V w (v = 0) + g 2 / \2A(3 A ) = f3~ A (-0.44 + 0.04# 2 ). 
The value of the effective potential in the minimum at v = tt — e is given by 



V R (v 



TT 



/T 4 (l8.0 - 0.26# : 



(6.4) 



for large g 2 (i.e. for e = 0.052). Thus we conclude that the minimum at v 7^ 
is deeper than that at v — for sufficiently strong coupling: g 2 > g 2 with the 
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critical coupling a c = g 2 /(4ir) « 4.9. Using the running coupling constant 
g 2 = -G2/ ln(A|fl 2 ) with the SU(N = 2) scale parameter A 2 and G 2 = 
47r(lliV/3) _1 , the following estimate is obtained for the critical temperature: 
T c = A 2 exp(G 2 /s c 2 )^A 2 . 

It is rather interesting that the non-trivial minimum of the effective po- 
tential due to the polynomic part Pq(v) is very close to v — ir corresponding 
to vanishing Polyakov line W = cos(t>/2) || |]. As discussed above, the ba- 
sic subprocess leading to this non-perturbative minimum is the decay of the 
A 03 condensate in two transverse gluons and a timelike fluctuation 5<p, which 
annihilate once again into the condensate. This subprocess is suppressed for 
v — » due to the infinitely increasing mass M ~ (1 — cosf ) -1 / 2 — > oo of the 
fluctuations S(f>, but it turns out to decrease the free energy of the vacuum for 
finite values of M, i.e. for v ^ 0. The decrease of the free energy is expected 
to have a maximum for the smallest value of the mass M, i.e. for v = 7r. In 
our result the small shift e of the minimum of the effective potential is most 
probably the consequence of the various approximations. 

7 Summary 

In the present article we have derived the renormalized finite temperature 
effective potential of the continuum SU(2) Yang-Mills theory as function of 
the non- vanishing expectation value v of the vector potential component A 03 , 
v = g(/3/a)(A 03 ). The quantum fluctuations around this mean field are taken 
into account up to the order g 2 perturbatively. The Haar measure has been 
taken into account in the definition of the partition function of the theory 
in order to treat large fields correctly. Due to the non-vanishing background 
field IR divergent terms occurred in the effective potential accompanied by 
UV divergences. Both are removed by choosing the IR momentum cut-off as a 
power series of the inverse of the UV momentum cut-off. The renormalization 
conditions requiring the cancellation of UV divergences have been solved in 
a linearized form. (Logarithmically divergent terms were neglected.) In the 
continuum limit the IR momentum cut-off turned out to be proportional 
to the temperature and takes values smaller than the first non-vanishing 
Matsubara-mode. 

In our approximation the main qualitative features of the finite tempera- 
ture effective potential are independent of the Haar measure. It reveals itself 
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only in the rather fine details. The renormalized effective potential has a min- 
imum at non- vanishing background field very close to v = it. For sufficiently 
strong coupling g 2 > g 2 , i.e. below the critical temperature T < T c « A 2 
(with the SU(2) scale parameter A 2 ) this non-trivial minimum is deeper than 
the minimum at v — describing the perturbative phase of the theory. It 
is well-known that a minimum at v — tt corresponds to vanishing Polyakov 
line, i.e. confinement. In our approach some important non-perturbative ef- 
fects are already taken into account by the mean field, although our treating 
the field fluctuations around the mean field perturbatively is not really jus- 
tified for strong coupling, Consequently, we cannot expect to obtain precise 
information on the details of the phase transition (the position of the non- 
perturbative minimum, the order of the phase transition, etc.) Nevertheless, 
we can conclude that the main mechanism responsible for the existence of the 
non-perturbative phase is the decay of the A 03 condensate (the background 
field v) in two transverse gluons and a timelike one. 

In order to get more reliable information on the position of the non- 
trivial minimum our treatment had to be improved at least by using the 
global center invariant tree level action proposed in [[H], |ll|], not to break 
global center symmetry explicitly at the beginning. Then all polynomials of 
v would be replaced by some polynomials of periodic functions of v and the 
non-trivial extremum would occur exactly at v = 7r. 
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A Contribution f\ 

The contribution fx is given by the expression (|3.1|). Using Eqs. ( p. 27] ) and 
( |2.30|) for the propagator we see that it is diagonal in the spatial indices (ij) 



for x = x', D l j£ b = ^8^D ab , and, consequently, the matrix D l j£ b = ^8 % ^D ab 
defined by Eq. ([3.2|), as well. Here we introduced the traces with respect of 
the spatial indices: 

D ab = jY,J-^[ Aab (^^) + 2A ab (u n ,p)] J (A.l) 
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and 



D 



ab 



a 2 . ^ 
7^ 



d 3 p 



A afe K,/i) + 2A a6 K,p) 



ab i 



with 



(A.2) 



(A.3) 



Performing the sums over the colour indices we obtain the following ex- 
pression: 



h = ^^D ll + D ll D^ + D ll D ll -D 12 D 12 ), (A.4) 



3a 4 

where the last term vanishes identically due to D 12 = 0. 
Further we make use of the following Matsubara sums: 

*(P»°) = ^£ rf n($ = ^(2n p + l), 



(A.5) 



*(p. ) = ^E^(P") = f (2^) + l) (A.6) 



with the Bose-Einstein occupation number n p = (e^ — l) and the 'gener- 
alized occupation number' N p (v) defined by 

^ «r / s , sinh(nV) . ( . 

2N p u +1 = wo ; ^ • A.7 

cosh(pp) — cosf 

Let us introduce the integral: 

R(v) = a 2 f j^v(p,v). (A ' 8) 

For v = we find: 
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with a = aA = (6tt 2 ) 1/3 - For v ^ let us write R(v) = R(0) + AR and 
<r(/i, v) = 0) + Act with 



AR = a 2 



~ a 



d 3 p sinh(/3p) 
(2tt) 3 2^ 
d 3 p /3p 



1 



1 



cosh(/3p) — cosf cosh(/3p) — 1 
1 1 



(2tt) 3 2p [ 1 - cos v + P V/2 /3 V/2 



27T/3 2 



sin - 
2 



and 



Aa = -(N^v)-n^). 



For the renormalization we need the limits 



i?oo(0) = lim R(0) = a 2 [ 



d 3 p 1 a 2 
(2itf2~p = 8^' 



and 



a 1 a oo (/i,0) = lima V(«, 0) = — . 

/3-*oo 2//a 

Furthermore we find A^i? = 0, i.e. Roo(v) = -Roo(O) and A^cr = 0. 
Putting the above formulas together we obtain for v = 0: 



(A.10) 



(A.ll) 



(A.12) 



(A.13) 



/i° 



v/ 



d 3 pd 3 q n p n q 



(2tt) 6 pg 
d 3 pd 3 q rip 
(2tt) 6 pg 



d 3 pd 3 q n p 
(2tt) 6 pg 



and 



36/3 4 



h 1 / 



(A.14) 



d 3 pd 3 q 1 
(2tt) 6 pg 
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+-(a(fi,0)-a oo (fi,0))a J j^yj 

2 r d 3 p 1 

+- (<r(/i, 0) - (Too^, 0) - a(ii, 0)) a 2 / 77—^7- 

+^a(fM, 0) {a(p,0)-2<T oo (p,0)) 



(A.15) 



leading to Eq. (|3.4|) . The terms linear in the occupation numbers cancel 
after renormalization in the glue gas contribution f® gl . 
The additional terms due to v 7^ are given by: 



Jl 3a 4 



4 ( -a(/i, 0) + 2R{Q) \ (-Aa + 2Ar\ 



-4 (-a^/i, 0) + 2R oo (0)) (-Aa + 2AR 
+ ^— Act + 2ARj 

-2-Aoocr (VV/i, 0) + 2R(Q) + -Aa + 2Ai?) 
a \ a J 



(A.16) 



Separating the contribution of the glue gas and that of the vacuum, we arrive 
at Eqs. (O) and (pi). 



B Contribution 



B.l Contracting the spatial indices 



By appropriate redefinition of the colour indices we obtain from Eq. (|3.14|) : 



3B 



abc gef 



-e e 



9^1 

■(27r) 3 (5(g + p + Q)f35 ni+n2+n3 ■ 
■a 2 \^D kieb (u ni ,q)D ljfc {uj n , 2 ,p) 

-D k ^ b (u ni ,q)D li ^(u n2 ,p}) Q i Q k D lj9a (u n3 , Q) 



^8 /33 



d 3 qd 3 pd 3 Q 



+ ~D 



kteb 1 
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+D lieb (u ni , q)D k ^ c (cu n2 ,p)) ptyD 1 ^^, Q) 

-D l ^ b (uj ni ,q)D k ^ c (uJ n2 ,p)) p i Q k D li9a (u n . v Q)] . (B.l) 

We start with contracting the spatial indices. For the sake of simplicity we 
write: U = A eb (oj ni , fi), U = A eb (uj ni ,q), V = Af c (uo n2 ^), V = A/ C (cu„ 2 ,$, 
W = Af a (cj„ 3 ,/i), W = A^(uj n . A ,Q). U, V, and W correspond to the longi- 
tudinal parts, U, V, and W to the transverse parts of the propagators. 

The expression 



~ q k q l 



hi 



(U-U)^ + U5 
q 2 

(U -U)^ + U5 kl 



(V -V)^- + V5 lj 
p 2 

l i 

(V -V)^+V8 li 

p2 



k I 

(U -U)(y -V)^(qY - qitf) 



+U(V - V)^(S ki p j - 5 fc V) + V(U - U)%r{5 l] q l - 5 li q j ) 
p 2 q 2 



+UV(5 ki 5 lj - 5 kj 5 li ) 
is multiplied with 



Q l Q k 



Q 



(w-w) 



(B.2) 



(B.3) 



in the first term of the bracket. The terms proportional to (W — W) vanish 
due to the contraction of a symmetric tensor with an antisymmetric one with 
indices ij, and we find: 



(U -U)(V -V)W 
+U(V-V)W[Q 2 - f 



(qQ) 2 (qQ)(qpWQ) 



q 2 p 2 



(PQ 



+ (U- U)VW2^Q^- 



-UVW2Q 2 . 



(B.4) 
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Proceeding similarly in the second and third terms we find for the bracket in 
Eq. O): 



uvwfaftftQ) + uvw<j) X {q,p,Q) + uvw<f> 2 (q,p,Q) + uvwMq,p,Q) 

(B.5) 



with 



Mq,p,Q) = Q -2{pQ) + — 2 ~ 



(qp)(qQ ) (qp)(pQ )(qQ ) , R „x 

+^ g 5-5 , (B.o) 



i 2 (9,P,Q) = ^^"2 ^ + ^2 > ( R? ) 



h{q,p,Q) = Q 5 5— + 5~2 ' ( R8 ) 



MV,P,Q) = — 2 ^2 • (B.9) 

The contribution f^s takes now the form: 

X 1 2 abc aef a2 1 V- [ d 3 qd 3 pd 3 Q 

/3B = "2^ e 6 ^ E 7 ( 27r )6 % + P + Q)^n 1+ n 2 +n 3 - 

A e Vm, q)Af c (u n2 ,p)A<> a (u n3 , Q)MQ,P, Q) 
+A eb K 1 , / i)A^ c K 2 ,^A^K 3 ,Q)0 2 (g,p,g) 
+A e6 (u; nu g)A /c K 2 , //)A^K 3 , Q) 
+A eb (u; ril ,^A fc (uj n2 , f i)A 9a (uj n3 ,Q)c/ )0 (Q,P,Q)] • (B.10) 
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B.2 Matsubara sums 

Let us now introduce the sum: 

S(q,P,Q) = 4 E ^„ 1+ n 2+ ns^e fle/ A c6 (w ni ,g)A^( £< ; n2 ,p)A» (a; n 3,^) 

P n\ri2nz 

E s B (q,p,Q) (B.ll) 

permutations of (qpQ) 

with 

Sb(q,P,Q) : 



- £ /55 ni+n2+ „ 3 rf- (p)rf„ 3 (Q). (B.12) 

r* ni?i2n3 

With the appropriate change of the momentum variables we can write: 
1 g 2 2 f d 3 qd 3 pd 3 Q 



hB — A 



— jO- 

2 a 4 



(2k) 

S(q,p,Q)(f) 3 (q,p,Q) 
+S (/i, p, Q)<f> 2 (q, P,Q)+S (q, fi, (g, p, Q) 

+5(/i,/i,Q)0 o (g,p,<5)] 

1# 2 2 f d 3 qd 3 pd 3 Q , 
~2^ a J ^2^ 6{q + P + Q) - 

■ s B (q,p,Q)P(f) 3 (q,p,Q) 

+ (s B (fi,p,Q) + s B (p,n,Q) + s B (Q,p,fJ.)) ®(q,p,Q) 
+2 (s B (/i, /i, Q) + s B (/i, Q, fi) + s B (Q, fi, /i)) <f> Q (q, p, Q) 



(B.13) 



with 



P03 



8(p 2 + Q 2 ) -2pg-6(pg) 



2 p 2 + Q 2 n (pQ) 

p 2 Q 2 



- 2- 



p 2 Q 2 



P 2 Q 2 ( P 2 + Q 2 ) + (p 4 + Q 4 )(pQ )-(p 2 + Q 2 )(pQ ? - 2(pQ y 



(p2 + Q2y _ A{ fQ )2 



(B.14) 
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(p 2 + Q 2 ) i- 



+2 



P 2 Q 2 
iir + Q 2 ) -1 



P 2 Q 2 



(p 2 + Q 2 )p 2 Q 2 + (p 4 + Q 4 )pQ - (p 2 + Q 2 )(pQ f - 2(pQ 
(p 2 + Q 2 ) 2 -A(pQ) 2 



(B.15) 



where we inserted g = — p — Q. 

We write for the Kronecker-delta: 



+n 2 +n 3 



e 0(g°+p°+Q°) _ 1 



(B.16) 



and perform the sum s# over the Matsubara frequencies by the contour 
integral technique of finite temperature field theory flUf : 



sb(q,p,Q) 



1 



h(q,p,Q) 



(B.17) 



with 



h(q,p,Q) 



n+ + l 



n- + 1 



+ 1 J(±£5=,±^,±^; 



(B.18) 



and i?^ = q±iv/{3. 

We express the sum h through the occupation numbers. Then we separate 
the vacuum term h , and the terms hi and hi linear and quadratic in the 
occupation numbers, respectively, h = ho + hi + hi with 



ho(q,P,Q) 



q + p + Q' 



(B.19) 



hi(q,P,Q) = ^riQ 



Q + P 2( + -x Q + Q 

(q + p) 2 -Q 2 {Up p> (q + Q) 2 -p 2 



+2{ < +n «\p + Q) 2 -q 2 - 



(B.20) 
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and 

h 2 ( q ,P,Q) = 4n Q « + V) g g4 - p4 - g4 q _ ~ q ^ p2 \ Q 2q2Q2 _ 2p2Q2 

p 2 _ g 2 _ Q2 

+4n Q (n g + n q )p ^ + p4 + Q4 _ 2q2p2 _ 2q2Q2 _ 2p2Q2 



+2 



71 



1 



n-)(n+ + n p )Q(Q 2 - q 2 - p 2 ) 



q 4 + p 4 + Q 4 - 2g v - 2g 2 g 2 - 2p 2 g 2 



(w+-7i g )(w+-7i p )ggp 
g 4 + p 4 + g 4 - 2q 2 p 2 - 2q 2 Q 2 - 2p 2 Q 2 ' 1 ' ' 

These expressions are symmetric in the momenta for vanishing background 
field v — 0, as expected. 



B.3 Subtraction of loop-integrals at zero temperature 

The 'renormalized' contribution f'3 is given in terms of the diagrams in Fig. 
4. As the momentum dependent functions 03, 0o, and $ are independent of 
the temperature, we can perform the subtractions on the sum s B directly. 
We have to subtract the limiting values of s B when two of the occupation 
numbers with all possible choices are taken in the limit (5 — > 00. 
If none of the momentum variables is equal to zero, we get: 

Wat® - i^+ vffi-^ W +">-)■ 

h<S?\q,P,Q) = ^\ n +2 p + ° ; «+",-) (B.22) 
and for the renormalized sum s(q,p, Q): 

s(q,p,Q) = s B ( q ,p,Q) + ^[h^\ q ,p,Q) + h^\q,p,Q) + ti 

n Q (n+ + n") q 2 - p 2 - Q 2 

P g g 4 + p 4 + g 4 - 2gV - 2g 2 g 2 - 2p 2 g 2 

K) + ^q)(< + ?V) q 2 -p 2 - Q 2 

ApQ q 4 +p 4 + Q 4 - 2q 2 p 2 - 2q 2 Q 2 - 2p 2 Q 2 
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(n+ - n p )(n^ - n Q ) Qp 

2Qp q 4 +p A + Q A - 2q 2 p 2 - 2g 2 Q 2 - 2p 2 Q 2 

' (B.23) 



2qpQ(q + p + Q) 



(The ~ sign means that s(q,p,Q) can be replaced in the integrand by the 
r.h.s. of Eq. ( |B.23|) .) Inserting q = —p — Q we obtain: 

1 pQ 
s(q,P,Q) = -^z(4n Q + n+ + n Q )(n+ + n p )- 



1 i 



8 V P W (PQ) 2 -p 2 Q 2 

1 

(B.24) 



2qpQ(q + p + Q)' 

The terms linear in the occupation numbers cancel due to renormalization, 
similarly as in the perturbative calculation without background field ||. 

If one or two of the momentum variables take the value of the IR cut-off 
fi, we make an expansion in its powers. The UV renormalization is performed 
by subtracting the corresponding limiting values with (3 — > oo for fixed finite 
fi. It removes the terms linear in the occupation numbers. Thus we obtain: 

s(p, fi, Q) + s(ji, p, Q) + s(Q, p, fi) = 

(4n Q + nJ +ng)(n+ +n~) p 2 + Q 2 



AQp (p 2 - Q 2 ) 2 

nQ(n+ + n-) n M (n J + Uq) 



Q(p 2 - Q 2 )/i 2Q(p 2 - Q 2 )/i 
ry(n+ +71") + n-){n+ + n~) 

2p(p 2 -Q 2 )fi 2pfi{p 2 -Q 2 ) 
2(n+ - n~) + n+ - n Q ] (n+ - n~) 



+ 



2(p 2 - Q 



2\2 



2pQ(p + Q) \fi p + Q / 



(B.25) 



s(fi, fi, p) + s{fi, p, fi) + s(p, fi, fi) = 
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2p 3 fi 
1 



(2n p + n+ + n p ) (nt + n u ) + (nj + n_ )n A 



Ap 2 fi 2 M 
1 



W + ™„ +4n M )(nJ + n ) 



2p 4 
3 



2p 2 /i 2 



1 - 



2/i X 



The 'renormalized' contribution / 3 can now be rewritten as 

g 2 fd 3 pd 3 Q_ T 



^ 2 J (2nY 

with the integrand 

2 = -s(q,p,Q)P(j) 3 (q,p,Q) 



(B.26) 



(B.27) 



- (s((i,p,Q) + s(p,p,Q) + s(Q,p,n))$(q,p,Q) 
-2 (s(fi, [i, Q) + s(fi, Q, p.) + s(Q, fi, fi)) <t> (q, p, Q) 



(B.28) 



and q = —p — Q. This can be cast in the terms of the order n°, n M , and 
n 2 ^ depending on the occupation numbers, and the term independent of the 
occupation numbers according to Eq. (|3.16|) . 

B.4 Momentum integrals 

Here is appropriate to make some remarks on the calculation of the momen- 
tum integrals. 

The momentum dependent functions multiplying the Matsubara sums 
in the integrands ( |3.17| )-( |3~2"3"D can be replaced by more simple expressions. 
They are obtained by separating the terms of the integrands with odd powers 
of the scalar product (pQ ) and leaving them out because they do not con- 
tribute to the integrals. Equations valid only up to those terms are denoted 
by ~. Thus the following expressions are found: 



p 2 Q 2 - (pQ Y p-Q 



2Q 2 



+ 



{P 2 -Q 2 Y{P A -Q A ) 



8Q 2 (p 2 + Q 2 ) 2 - A(pQ ) 2 



(B.29) 
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(PQ)pMq,p,Q) 



(pQ) 2 -p 2 Q 2 



, (pQ) 2 

' p 2 Q 2 



4 + 



(p 2 + Q 



2\2 



p 2 Q 2 



4 + 



{p 2 + Q 
P 2 Q 2 



2\2 



(p 2 + Q 



2\2 



(p2 + Q2)2 _ 4(^g ): 



(B.30) 



(pQ ) 2 - p 2 Q 2 



P 2 + Q 2 
P 2 Q 2 



p4 + g4 + 6p 2 Q 2 

2( P 2 + Q 2 )pQ n 



(B.31) 



(obtained after performing the angle integration) , 

2$(g,P,Q) 



2(p 2 + Q 2 )(p 2 -Q 2 ){-^ 



Q 2 



+ 



1 + 



(p 2 + q 

4p 2 Q 2 



2\2 



(p2 + Q2)2 _ 4( ^Q )2 



(B.32) 



There is a remarkable cancellation of several terms in the integrand of X 2 
^see Eq. ( frTTD ): 



(p 2 - Q 2 ) 2 



~ 5 2(/ + Q 2 ) 2 



2 {pQ) 2 -p 2 Q 2 
(P 2 + Q 



2\2 



+ 1 + 



($3 ) 2 

+2^+ + 2 | 1 + 



4p 2 Q 2 v 
(p 2 + Q 2 ) 2 



4p 2 Q 2 J (p2 + Q2)2 _ 4 (^Q )2 



p 2 Q 2 



4p 2 Q 2 



-2 1 + 



Cp 2 + Q 

o2/02 



2\2 



(p 2 + Q 



2\2 



4p 2 Q 2 J (p 2 + Q 2 ) 2 - 4(pQ )2 



, M (pQ ) 

p 2 Q 2 3 



(B.33) 



The last equation is obtained by performing the angle integration. 

For evaluating the glue gas contribution the following approximation 
is used: 



Up +n p - 2n p 



+ 



I3p p -iv 
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2 (e^cosu - 1 
e 2i3 P _ 2 e Pv cos v + 1 - 1 

f for p> 1/(3 
\ v{jp) for p < 1/(3 



where 



cos v — e P p 



1 - 2e~PP cos v + e- 2f3 P 1 - 
Then we obtain from Eqs. (ft. 17 ), and (ft.24 ): 



^B.34) 



(B.35) 



& « 2 7^)4 jf ^ / M [^(p) + "(pMQ)] • (B.36) 

We use the saddle point approximation and replace the expression pv{p) by 
its value taken at po = 1/(3: 



G(v) = (3p Q v(p G 



cosf — (1/e) 



1 -2(l/e) cost; + (1/e 2 ) 1 - (1/e) 



(B.37) 



This leads to the expression (ft.28|) for f£,. 
For the evaluation of the integral 



y,2 



^1 



g d pd Q-(n+ -n p )(n^-n Q ) 



{pQ) 2 -p 2 Q 2 



(B.38) 



we make use of the approximation: 

1 1 



n p -n p 



e f3p+iv _ e /3p-iv _ I 



for p > 1/(3 
5(p) for p <l/(3 



(B.39) 



with 



-/3p 



-0p p—iv 



-Pp 



—l smv 



^l/e)(cosh 1 — cosv) 

(B.40) 



Then we obtain the expression ( ft.29| ) for /g 2 . 
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The contribution 

fl R = fsi + U + M + fii + fap (B.41) 

was estimated in the following way: (i) the integrals containing occupation 
numbers in their integrands were calculated by using the low-momentum 
(p, < Q < 1/(3) expansions 

n Q « — - - + — PQ - —WQf + WQf + ..., (B.42) 

Q PQ 2 12 PV 72(T v; 720-42 VPV; ' v ; 



n+ + riQ « _i + A/?g + ^(/3Q) 3 + T ^C(/3Q) 5 + ..., (B.43) 



» - i Asmv(l-±A(f3Q) 2 + ±D(f3Q) i + ..)) (B.44) 

with A = (1 - cos v)-\ B = A - 3A 2 , C = A - 15A 2 + 30A 3 , D = —A + 6A 2 ; 
(ii) for large momenta (1//3 < Q < A) the estimates nq « e _/3< ^, + Hq « 
_ e -pQ ) anc [ n + _ ~ —1^4 sin f e _/3< ^ are used; fnij slowly varying factors 
in the integrands were generally separated and taken out from the integral 
with their value at the maximum of the rapidly varying factor; (iv) all IR 
divergent terms were included. 



C Contribution j L 



4 



Calculating the bare contribution f^s we performed the sums over colour 
indeces at first. The product of propagators was rewritten as 

D Hab (u n ,p)D Ma ' b '(uj m ,q) = 
= A ab (uj n ,p)A a ' b '(uj m ,p)cos 2 6 



+ 



A ab (iu n ,p)A a b (u m , p) + A ab (u n , p) A a b (u m , q) 



sm 2 9 



+A ab (uj n ,p)A 



a'b' 



cos 2 fl 



(G.1) 



where 9 is the angle of the vectors p and q. Making use of expression ( 2.28ft 
for the matrices A ab (uj n ,p) and the fact that the sums contain only the terms 
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with 112 = —TI3, and d^ip) = d^ n (p) we obtained: 



' a 4 (3 2 



4B 



U2+U3 



n2«3 



d 3 p2d 3 pr 



V(~P2 - Ps) 



(2tt) 6 

■{(u n3 -(v/f3)) 2 d± n2 (0)d+ 3 (0)cos 2 e 
+2K 3 -W/3)) 2 dl n2 (p 2 )< 3 (0)sin^ 
+ K 3 - (t;//?)) 2 dl n2 (p 2 )<(p 3 )(l + cos 2 (9) 
+2(^//3)a; n 3< (p 2 )< (p 3 ) (1 + cos 2 0) } 



= -2 



^7 



d 3 pd 3 q 



(2tt) 



V(-p-q) 



■ s 2 b(v,v) cos 2 6 + 2s 2 b(p, A 4 ) sin 2 6* 
+s2b(p, g)(l + cos 2 0) + 2si B (p, + cos 2 0) 
with the Matsubara sums: 



(C.2) 



SiB(p,q) 



E 



/? 2 [(P° + i(v/[3)f - P] [(^ + i^//?)) 2 - p 
1 



% < + l 1 n+ + 



h{±E;,±Ej), 



(C.3) 



and 



S2B(p,q) = 



1 _ e /3(p°+?°) _ 1 



_ y 

/9 2 4^3 P° + ?° + (^//?)) 2 + p 2 + (v/(3)f + q 2 



V - («/0)) 5 



= --T 



h(p°,q°) 



P 2 ntT 3 [(p° - i(^//5)) 2 - P] [(q° + i(v/[3)f - q 2 



where 



e f3{p°+q°) _ 1 



p° + q 



— (-i(^)gV), 



(C.4) 



(C.5) 
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WV) = ' p0+ ~ > " -K»/fl) V, (C6) 

and n± = (e^* — l) , -E^ = p±i(v/@). Let us express these sums in terms 
of the occupation numbers: 

4pgsi B (p,g) = (- - , „ 1 o; ^, + — ] — )iCE-a\ 



p + q-2iC -p + q-2iC 



+ \p-q-2iC + p + q 1 +2ic) 
~ { p + q-XC + p-q-Tic ) 
+ { p-q + TiC + p + q l +2ic ) 
+1 ^ 2 (^y-^^)' (C - 7) 



4pqs 2B (p,q) = ^ P , \(q 2 - 4C 2 )(n+ + n q ) + i4gC(n+ - n g ) 



^2 _ g2 

-^\(? 2 - 4C 2 )« + V) - -|^i4gC(< - V) 
" 2 -(q 2 -AC 2 ). (C.8) 



p + q 

The subtractions shown in Fig. 5 can be carried out on the sums Slb(p, q) 
and S2b(p,q) directly. Then the 'renormalized' contribution / 4 must be cal- 
culated by replacing the bare sums via 

Sa = S aB - S M- S M- 8 W (C.9) 

(a = 1,2), where the upper indices (p), (g), and (jog) denote that the occu- 
pation numbers , n^, and both and are taken in the limit f3 — > oo, 
respectively. As a result the terms linear in the occupation numbers cancel 
in the 'renormalized' expressions and we obtain: 

ipM = WIV ^^^, , (CIO) 
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4pqs 2 (p,q) 



-2- 



p + q 



(q 2 - 4(v/[3f 



(C.ll) 



Replacing the bare expressions by these 'renormalized' ones we obtain: 



2g 2 f d 3 pd 3 q 



—a 



— eT 



,a/x 
p 2 



(27T) 
/ 

1 



M 2 + (p + g) 2 
'a//3) 2N 



4a/i 



4tT 



(a/i) 2 



cos 2 fl 



(a/i) 2 , 
(V - 4( W //3) 2 ) (l + cos 2 fl) 



2pg(p + g) 

P-g (^) 2 
pq (p + q) 2 + 4(v / (3)' 



(l + cos 2 6 



(C.12) 
oo. Therefore 



For v — the contribution (|C.12j ) vanishes due to M 
this is a vacuum contribution for v ^ 0, — f± vac - 

The momentum integrals are treated similarly as those in Appendix B. 
The integral over the angle 9 is performed at first. Then estimates are ob- 
tained by casting the integrands in slowly and rapidly varying factors and 
taking the slowly varying factor out of the integral with its value at the max- 
imum of the rapidly varying one. Keeping the terms up to the order g 2 and 
making use of M 2 ~ g 2 A 3 , we obtain the final expression (|3.32|) . 
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